Narrow-Band Biphoton Generation near Atomic Resonance 
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Generating nonclassical light offers a benchmark tool for the fundamental research and potential 
applications in quantum optics. Conventionally, it has become a standard technique to produce the 
nonclassical light through the nonlinear optical processes occurring in nonlinear crystals. In this re- 
view we describe using cold atomic- gas media to generate such nonclassical light, especially focusing 
on narrow-band biphoton generation. Compared with the standard procedure, the new biphoton 
source has such properties as long coherence time, long coherence length, high spectral brightness, 
and high conversion efficiency. In this paper we concentrate on the theoretical aspect of the en- 
tangled two-photon state produced from the four-wave mixing in a multilevel atomic ensemble. We 
show that both linear and nonlinear optical responses to the generated fields play an important role 
in determining the biphoton waveform and, consequently on the two-photon temporal correlation. 
There are two characteristic regimes determined by whether the linear or nonlinear coherence time 
is dominant. In addition, our model provides a clear physical picture that brings insight into under- 
standing biphoton optics with this new source. We apply our model to recent work on generating 
narrow-band (and even subnatural lincwidth) paired photons using the technique of electromag- 
netically induced transparency and slow-light effect in cold atoms, and find good agreements with 
experimental results. 
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1. Introduction 

Nonclassical light generation has attracted much atten- 
tion over the last 40 years, partly because it not only 
provides a powerful probe for addressing fundamental is- 
sues of quantum theory such as complementarity, hid- 
den variables, and other aspects central to the founda- 
tions of quantum mechanics [1]; but also holds promise 
for many potential applications to quantum information 
processing [2], quantum computation and communica- 
tion [3], quantum cryptography [4], quantum imaging 
[5], quantum lithography [6, 7], and quantum metrol- 
ogy [8]. In particular, entangled photon pairs have al- 
ready been established as a standard research tool in the 
field of quantum optics. Traditionally, paired photons 
are produced from spontaneous parametric down conver- 
sion (SPDC) [9, 10], in which a strong pump laser drives 
the atomic oscillators in a noncentrosymmetric crystal 
into a nonlinear regime, and then two down-converted 
beams are radiated by these oscillators. The two down- 
converted photons from such a nonlinear process usually 
have very broad bandwidth (typically in the THz range) 
and very short coherence time (typically around few pi- 
coseconds) so that their waveforms are not resolvable by 
existing single-photon detectors (which have resolution 
around the nanosecond range). In addition, the short 
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coherence length of the SPDC biphotons, roughly about 
~ 100 /jm, also limits their application in long-distance 
quantum communication. The optical properties of these 
wide-band biphotons have been well studied both exper- 
imentally and theoretically for more than two decades 
[11, 12, 13]. 

Narrow-bandwidth biphotons are ideal for a number of 
recently proposed protocols for long-distance quantum 
communication based on coherent interaction between 
single photons and atomic ensembles [14, 15], which re- 
quire efficiently absorbing biphotons and storing the en- 
tanglement. Other applications include the long-distance 
quantum state teleportation which requires long tempo- 
ral coherence time of paired photons [16]. Using cavity- 
enhanced SPDC in nonlinear crystals, paired photons 
with a bandwidth of about 10 MHz and coherence time of 
about 50 ns have been reported [17, 18, 19, 20]. However, 
the SPDC process is usually tuned to occur in the region 
of far-off- resonant atomic transitions, which limits the 
conversion efficiency and requires a large pump power. 
Moreover, in the far-off-resonant pumping the nonlinear 
parametric coupling coefficient can be treated as a con- 
stant and the information about the material internal 
transition structure is unresolvable from the two-photon 
correlation measurement. 

The two most recently developed advanced technolo- 
gies make the near- and on-resonance nonlinear processes 
possible. The first is laser cooling and trapping of neu- 
tral atoms [24]. Because cold atoms with a temperature 
of about 100 fi K have negligible Doppler broadening and 



very long lifetime, their atomic hyperfine structures can 
be resolved without need of a Doppler-free setup. The 
second is the electromagnetically induced transparency 
(EIT) [21, 22, 23] which not only can eliminate absorp- 
tion on resonance but also enhance nonlinear interactions 
at low light level [25, 26, 27, 28]. Following the theoretical 
proposal in [14] and early experiments by Lukin and Kim- 
ble groups [29, 30], generation of nonclassical correlated 
photon pairs by exploring the "writing-reading" proce- 
dure has attracted much research interest. However, the 
created paired photons are not maximally time- frequency 
entangled due to two separate operation processes. Using 
cold atomic-gas media and EIT-assisted four-wave mix- 
ing (FWM), Harris group at Stanford University was the 
first to generate time-frequency entangled narrow-band 
biphotons with cw driving fields in a double-A system 
[31, 32]. By controlling the EIT at a sufficiently high 
optical depth (OD) around 50, biphotons with subnatu- 
ral lincwidth have been produced in a recent experiment, 
where the coherence time is up to about 1 /is in a two- 
dimensional magneto-optical trap (MOT) [33]. Vuletic 
and his colleagues have also produced biphotons with a 
temporal coherence length of 100 ns by trapping cold 
atoms inside an optical cavity [34]. Du et at, by driving 
a two-level system into the nonlinear region, have also 
created paired photons with correlation time of same or- 
der of magnitude [35] . 

To describe the physics behind narrow-band biphoton 
generation from the EIT-based multilevel atomic system, 
two different approaches have been presented in the lit- 
erature. One approach is to use the Heisenberg-Langevin 
method by solving the coupled field operator equations 
in the Heisenberg picture, and the results agree well with 
the experimental data in the two-photon correlation mea- 
surement [36, 37]. The advantage of this method is that 
it provides a sophisticated calculation of the effects from 
the Langevin noises, the medium absorption and optical 
gain on the biphoton wavepacket. Working in the inter- 
action picture, the other approach developed by Wen et 
al. [38, 39, 40, 41] focuses on the description of the two- 
photon state vector after the nonlinear medium using 
perturbation theory. This state vector theory not only 
properly describes the experimental results but also offers 
a clear physical picture on biphoton generation mecha- 
nism. It is the intent of this paper to review and gen- 
eralize this latter approach by taking into account both 
linear loss and gain, and provide more insights on the 
roles of nonlinear interaction and linear response to the 
biphoton amplitude. Using this state vector picture, we 
show that the two-photon wave function is a convolution 
of the nonlinear and linear optical responses. According 
to this convolution, the two-photon temporal correlation 
is considered in two regimes, damped Rabi oscillation and 
group delay. 

The paper is organized as the following. In Sec. 2, us- 
ing perturbation theory we lay out a general formulism of 
the two-photon state and biphoton wavepacket (or am- 
plitude) generated from the FWM parametric process. 
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Fig. 1. (color online) Biphoton generation via a four- level 
double-A atomic system, (a) The level structure, where 
in the presence of a cw pump (ui p ) and coupling (ui c ) 
beams, paired Stokes (ui s ) and anti-Stokes (uj as ) photons 
are spontaneously created from the four-wave mixing pro- 
cesses in the low-gain regime, (b) The forward generation 
configuration, (c) The backward generation geometry. 



In Sec. 3, taking a four- level double-A atomic system 
as an example, we show that the linear and nonlinear 
optical responses to the generated fields play an impor- 
tant role in determining the two-photon amplitude. By 
looking at the nonlinear susceptibility, we illustrate the 
mechanism of biphoton generation in such a system. In 
Sec. 4 we describe the damped Rabi oscillation regime 
where the biphoton amplitude is mainly characterized 
by the structure of the third-order nonlinear susceptibil- 
ity. We show that the two-photon temporal correlation 
can exhibit damped and over-damped Rabi oscillaitons. 
In Sec. 5, we describe the group-delay regime where the 
biphoton wavepacket is mainly determined by the phase 
matching due to the linear optical response. In Sec. 6, 
we examine the two-photon interference using biphotons 
generated from such a cold atomic-gas medium. We show 
that by manipulating the linear optical response, it is 
possible to switch the two-photon anti-bunching-like ef- 
fect to bunching- like effect. Finally, we draw the conclu- 
sions and outlook in Sec. 7. 

2. The Two-Photon State and Wave Function 

A schematic of biphoton generation via a four-level 
double-A atomic system is shown in Fig. 1, where in the 
presence of a cw pump (u> p ) and coupling (w c ) lasers, 
phase-matched, paired Stokes (uj s ) and anti-Stokes (uj as ) 
photons are spontaneously produced from the FWM pro- 
cess in the low-gain limit. As shown in Fig. 1(a), the 
strong coupling laser forms a standard three-level A EIT 
scheme with the generated anti-Stokes field. Therefore, 
the role of the coupling laser here is that it not only 
assists the FWM nonlinear process but also creates a 
transparency window for the anti-Stokes photons by the 
slow-light effect. With the use of EIT, it is possible to 
generate the anti-Stokes photons near atomic resonance 
|1) — > 1 3), which greatly enhances the efficiency of FWM 
process. The phase-matching condition allows forward 



and backward two generation configurations, as shown in 
Fig. 1(b) and (c), respectively. In the backward genera- 
tion configuration paired Stokes and anti-Stokes photons 
can propagate collincarly with the pump and coupling 
beams [31, 33], and can also propagate in a right-angle 
geometry [32]. These flexible generation setups depend 



on whether the relationship k p 



is well satis- 



fied, where k PyC are wave vectors of the input pump and 
coupling fields. One advantage for the backward genera- 
tion geometry is that it allows us to easily separate the 
generated weak fields from the strong inputs. 

In the following discussions, we assume that the cold 
atomic-gas medium consists of identical four-level atoms 
prepared in their ground state |1) [see Fig. 1(a)]. The ide- 
alized atoms or molecules are confined within a long, thin 
cylindrical volume with a length L and atomic density is 
N . For simplicity, in this paper we will not take into 
account the Doppler broadening and polarization effects. 
In the two-photon limit, the quantum Langevin noise in- 
troduces unpaired photons which arc not of interest here 
and so is ignored. In addition, we concentrate on the 
two-photon temporal correlation. 

We start with the geometries shown in Fig. 1(b) 
and (c), where Stokes-anti-Stokes photon pairs propa- 
gate along the z axis. We denote the electric field as 
E = ±[£ (+) + £ (-) ] = \[E (+) + c.c], where #(+) and 
E^> are the positive- and negative-frequency parts. The 
input pump and coupling beams are assumed to be strong 
classical fields, so that 



£^ +) (z,£) = E p e i[±kpZ - UJ " t] , 
E { +\ Zl t) = E c e i[kcZ ~^ t] , 



(1) 



where k PyC are pump and coupling field wave numbers. 
The ± sign in Eq. (1) stands for the forward and back- 
ward propagation geometries for the pump laser. The 
single-transverse-mode Stokes and anti-Stokes fields are 
taken as quantized, 
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where k s ^ as are wave numbers of Stokes and anti-Stokes 
photons, and A is the single- mode cross section area.. 
The annihilation operators a s (w s ) and a as (uj as ) satisfy 
the commutation relation 

[o,(«), oJ(c/)] = [&„.(«), &L(w')] = 6(u - «')• (3) 

In the interaction picture, the effective interaction Hamil- 
tonian for the four-wave mixing parametric process takes 
the form [38, 39, 40, 41] 



Hi = «£ [ L/2 ^^^H + h.c, (4) 

4 J-L/2 



where h.c. means the Hermitian conjugate, and \ IS 
the third-order nonlinear susceptibility to the Stokes (or 
anti-Stokes) field defined by the nonlinear polarizability 



P t 



(3)(+) 



eoX {3) E ( p +) Ei +) ELl [42]. With use of Eqs. (1) 
and (2), after the z integration we rewrite the Hamilto- 
nian (4) as 

Hi = — — I du as duj s K(u as , Lj s )sinc(——) 

x al(oo a8 )at(w s )e- i ^ + ^-^-^ t + h.c, (5) 

where Ak = k as ± k s — (k c ± k p ) is the phase mismatch- 
ing for the forward (+) and backward (— ) configurations. 
When the pump and coupling fields are incident with an 
angle respect to the z axis, the phase mismatching can 
be generalized as Afc = (fc as + k s — k c — k p ) ■ z, where 
z is the unit vector along the +z axis. The integral 
over the length L gives a sine function sinc(^|^), which 
determines the two-photon natural spectrum width. 
n{uj as ,uj s ) = -i ^™™™ s ^(3) (u) a s,Us)E p E c is the nonlin- 
ear parametric coupling coefficient, where zu as and vo s are 
central frequencies of the anti-Stokes and Stokes fields, 
respectively. 

Perturbation theory gives the photon state at the out- 
put surface(s) approximately as a linear superposition of 
|0) + I*), where |0) is the vacuum state. The two-photon 
(biphoton) state |\I/) can be expressed as 



I*) 



+ DC 



dtH^O). 



(6) 



Because the vacuum is not detectable, from now on we 
will ignore it and only consider the two-photon part. 
Higher order terms with four, six, etc. photons are neg- 
ligible for low-power continuous pumping. Using the 
Hamiltonian (5) in Eq. (6), the time integral yields a 
5 function, 2ttS(lo c + lo p — u> as ~ ^s), which expresses the 
energy conservation for the process. If L is infinite, then 
the sine function Lsinc(-^|^) in (5) becomes a 8 function, 
<5(Afc). In this case, the conditions 



LU C + LU p - LU as - 0J S = 0, 

kas ^ k 8 \k c i k p ) = u 



(7) 



both hold and the phase matching is said to be perfect. 
The phase matching condition arises from the fact that 
the FWM process is a coherent process. Now the two- 
photon state (6) becomes 

I*) = L I du as K(uj as ,u; p + lo c - uj as )s\nc{— — J 

xal s (uJa S )al(ujp +u c - w o «)|0). (8) 

As seen from Eq. (8), the two-photon state is entangled 
in frequency and wave number, but is not entangled in 
polarization. In frequency space, the entanglement is the 
result of the frequency phase-matching condition, which 
implies that the detection of a photon at frequency uj as 



requires the detection of the other photon at frequency 
Lu p + uj c — Lu as . The frequency correlation has interest- 
ing consequences for the temporal behavior of the pair 
[31, 32, 33, 36, 37, 38, 39, 40, 41, 44], as we shall sec. The 
state is also entangled with respect to the wave number 
since the sine function sinc(^|^) cannot be factorized as 
a function of k as times a function of k s . In the general 
noncollinear case, the wave-number entanglement has im- 
plication for the spatial correlation of photon pairs [38] . 
The polarization-entanglement generation has been dis- 
cussed in [40, 43]. 

To study the optical properties of paired Stokes and 
anti-Stokes photons generated from such a four-level sys- 
tem, let us look at the simple experiments of two-photon 
joint-detection measurement as shown in Fig. 1(b) and 
(c), where anti-Stokes photons go to detector Di and 
Stokes photons go to D2 . The two-photon temporal cor- 
relation is the quantity of primary interest in this paper. 
The field annihilation operators in time domain on the 
surfaces are expressed as 

a.(t) = -tL j dLoa s {Lo)e^ k ^ L ' 2 -^\ 



a as (t) = -^= [ dua a .(u;)e i V t "M L ' 2 - ut l (9) 

V 27T J 

It can be shown that they obey the following commuta- 
tion relation 

[a s (t),&l(t')] = [a a s(t),al(t')} = 6(t-t'). (10) 

The joint-detection probability is defined by Glauber's 
theory, 

G (2) (t as ,t s ) = (y\al(t s )al s (t as )a as {t as )a 8 (t s )\y) 
= \(Q\a as (t as )a s (t s )\V)\ 2 
- |*(W S )| 2 - (11) 

In general, ^(tas, t s ) is referred to as the two-photon am- 
plitude, or biphoton wavepacket. With the help of Eqs. (8) 
and (9), the two-photon wave function (or amplitude) 
now is 



*(i os ,i s )=V(r)e-^ + ^^, 



(12) 



where 



V>(r) = ^ J duj asK (u as )^(io as )e- l ^ T , (13) 

with the relative time delay r = t as — t s and n{uj as ) = 
n(uj as , lo c + Up — w a s)- In Eq. (13), <&{uj as ) is defined as 
the longitudinal detuning function, 



<b(u> as ) = sine 



/AkL\ 



,i(k as +k s )L/2 



V 2 / 



(14) 



From Eq. (13) we see that the biphoton wave function 
is determined by both the nonlinear coupling coefficient 



(k) and the longitudinal detuning function ($). In the 
time domain, ip(r) is a convolution (*) of k and $ 

i/)(t) = L[k(t) * <£>(r)] = L J dT'k{T r )^{T - t'). (15) 



Here, k(t) is Fourier transform of the nonlinear coupling 
coefficient K(u) as ), 



k(t) = — / duj as K(u) as )e 

Z7T 



— tUJ as T 



(16) 



and similarly, $(r) is Fourier transform of the longitudi- 
nal detuning function $(w as ), 

Ht) = ^J dw as $(w as )e- iu "° r . (17) 

The photon pair emission rate R can be evaluated, by 
using Eqs. (13), as 

R= h I ' d ^M^asMuj as )\ 2 = f dr|V-(r)| 2 . (18) 

Now let us consider a photon counting experiment. 
Let the detectors be at equal distances from the source 
so that the relative time delay caused by the free-space 
propagation of electric fields are equal. Suppose the 
Stokes photon is detected by D2 at time t s , and after 
a controllable delay (r) detector Di is opened to detect 
the anti-Stokes photon. Assuming perfect detection ef- 
ficiency, the averaged two-photon coincidence counting 
rate is determined by 

Rcc{t) = lim — / dt s / dt as 
T-,00 1 J J 

X G^(t as ,t s )Il(t as -t s ;T,T + t c ) 

dr / G< 2 V)n(r';T,T + tc). ( 19 ) 



II(r'; t,t + t c ) is the coincidence window function within 
the time bin width t c : II = 1 for r < r' < r + t c ; oth- 
erwise, II = 0. When the two-photon coherence time is 
much larger than the detector time bin width t c , Eq. (19) 
becomes R cc (t) <~ G( 2 )(r)t c , which is a measure of the 
two-photon correlation function. This is the essence of 
narrow-band paired photons measurement discussed in 
this paper. 

Equation (12) is the starting point for the discussions 
in the following sections. By extending wave numbers as 
complex variables, the effects of loss and gain inside the 
material can be properly included in the above theory, 
as we shall see. It is found that to be consistent with 
the Hciscnberg-Langevin theory [36, 37] in the low-gain 
limit, the argument in $ should be replaced by Afc = 
(k as + k* — k c — k p ) ■ z, where k* is the conjugate of k s . 

Before proceeding the discussion let us make a short 
summary. In this section, we have laid out the general 
formulism of the two-photon state and wave function gen- 
erated in a FWM medium pumped by two cw classical 



fields using perturbation theory. The two-photon state 
(8) is a coherent superposition and it is entangled in fre- 
quency and wave numbers. The biphoton wave function, 
as shown in Eq. (15), is a convolution of R and $. Alter- 
natively, Eq. (13) states that the two-photon spectrum 
is determined by both the nonlinear susceptibility and 
the phase matching. This implies two regimes for the 
temporal correlation measurement, as discussed in Sec. 4 
and Sec. 5. The methodology explored here is equiva- 
lent to the Heisenberg-Langevin theory [36, 37] at the 
two-photon level. Moreover, this method offers a clear 
physical picture of the paired photons, which makes it 
easier to understand the two-photon interference experi- 
ment [45, 46, 47, 48]. Using the above theory new phe- 
nomena have been analyzed [33, 44, 49] which are not 
easy to be interpreted from the Heisenberg-Langevin the- 
ory For example, the two-photon correlation at a high 
optical depth [33], the simulation by numerically solv- 
ing the Heisenberg-Langevin coupled-operator equations 
shows a good agreement in the leading edge spike with 
the experiment data, but the nature of optical precursors 
is much easier to be understood in the present biphoton 
state picture [44]. 

3. Nonlinear and Linear Optical Responses 

As shown in Eq. (12), the two-photon wavefunction is 
determined by both the nonlinear susceptibility and the 
longitudinal detuning function. Therefore, in this section 
we look at the linear and nonlinear optical responses to 
the generated fields. On the one hand, the importance 
of the nonlinear susceptibility is that it not only charac- 
terizes the strength of the nonlinear parametric process 
but also is the mechanism of biphoton generation. As a 
consequence, it determines the feature of the two-photon 
temporal correlation in one regime. On the other hand, 
the longitudinal detuning function (related to the linear 
susceptibilities) sets the natural spectral width of bipho- 
tons, and governs the pattern of the temporal correlation 
in another regime. The optical response of a four-level 
doublc-A EIT system has been discussed in [36, 37, 38]. 
Wen et al. have further systematically studied the opti- 
cal responses in a three-level double-A scheme [40] and a 
two- level system [39] . Recently, the responses of a four- 
level inverted- Y system have been analyzed in [41]. 

For simplicity, here let us still choose paired Stokes- 
anti-Stokes generation shown in Fig. 1 as an example to 
illustrate the physics behind FWM. Assuming that the 
intensity of the off-resonant pump laser is chosen so that 
the atomic population remains primarily in the ground 
level |1), the third-order nonlinear susceptibility for the 
generated anti-Stokes field is [28, 38] 



X 



(?>) 



N Mi3M32M24M4i /(eo?i 3 ) 



(A p + i7i 4 )p c | 2 - 4(w + Z7i 3 )(w + ryi2)] 

-A^13M32M24M4l/(eofi 3 ) 

4(A P + i7i 4 )(w - fl e /2 + ij e )(u> + VL e /2 + ij e ) 



(20) 



are dephasing rates, respectively. A p = lu p — um is the 
pump detuning from the atomic transition |1) — > |4). 
lu = uj as — 0J31 is the detuning of the anti-Stokes pho- 
tons from the transition |1) — > |3), and we will take 
Was = W31 as the anti-Stokes central frequency. Q e = 
^/|fi c | 2 — (713 — 712) 2 is the effective coupling Rabi fre- 
quency. 7 e = (712 + 713) /2 is the effective dephasing 
rate. The third-order nonlinear susceptibility has two 
resonances separated by O e and each is associated with 
a linewidth of 27 e . The two-photon coherence time is 
inversely proportional to 27 e . More importantly, the two 
resonances here indicate two types of FWM processes. 
Alternatively, two types of paired Stokes and anti-Stokes 
photons can be generated from these two FWM processes 
in the spontaneous emission region. One FWM process 
happens when the central frequency of the anti-Stokes 
photons is zu as + Q e /2 and the central frequency of the 
correlated Stokes photons is W42 — A p — fi e /2. The other 
FWM occurs as the anti-Stokes photons have the central 
frequency at w as — f^e/2 while the paired Stokes photons 
have the frequency at u>42 + A p + f2 e /2. This is the gener- 
ation mechanism of biphotons produced from the FWM 
in atomic-gas media. Consequently, the interference be- 
tween these two types of biphotons will appear in the 
two-photon temporal correlation, as shall be discussed 
in Sec. 4. The results obtained here agree well with the 
dressed-state picture [39, 40, 41]. 

The linear susceptibilities at the Stokes and anti-Stokes 
frequencies are, respectively, 



X«M 



N\^24\ 2 (u;-i"f 13 )/(eoh) 



W 



\Q C \ 2 - 4(lu - »7i 3 )(w - 4712) A2 + 7l 2 4 



, (21) 



Xas{u) = 



4N\fi 13 \ 2 (Lo + ij 12 )/(eoh) 
IO:! 2 - 4(w + i7i3)(w -M712) 



(22) 



In Eq. (20), /i^ are the electric dipole matrix elements, 
f^c = l^2zE c /U is the coupling Rabi frequency, and 7^ 



where Q, p — /inEp/h is the pump Rabi frequency. The 
complex wave numbers of Stokes and anti-Stokes photons 
are obtained from the relations k s — (uj s /c)y/l + Xs & n d 
k as = (u> as /c)y/l + \as, where the imaginary parts stand 
for the Raman gain and EIT loss, respectively. The linear 
susceptibility (22) to the anti-Stokes field takes the form 
of standard EIT. As expected, the coupling beam £l c not 
only assists the propagation of the anti-Stokes photons 
by creating an EIT window, but also manipulates the 
group velocity using the slow- light effect. Another role 
played by EIT here is that it allows the nonlinear optics 
occurring near atomic resonance with small absorption 
and hence enhances the efficiency of the nonlinear inter- 
action. 

It is instructive to examine the linear responses un- 
der some simplifying approximation. Taking \Cl p \ -C A p , 
Eqs. (21) and (22) give Xs — and k as ~ k asa +cj/V g + ia 
so that the wave-number mismatching is Afc ~ u>/V g +ia. 
Here fc as0 is the central wave number of the anti-Stokes 
field, V g is its group velocity, and a, the imaginary part 
of the anti-Stokes wave number, characterizes the EIT 



loss, respectively. Now Eq. (14) can be approximated as 



( u>L aL\ ( luL aL 

+ z— J cxp U 



$(w) " SinC fe •'2J~~V2V t 



2 ), (23) 



which sets a coherence time about t„ 



L/V g 



In 



Eq. (23), a = 27Vcri 3 7i27i 3 /(|^ c | 2 + 4712713) where 
<7i3 = 27r|^,i3| 2 /(£:o?iAi37i3) is the on-resonance absorp- 
tion cross section in the transition |1) — > |3). 

As discussed in this section, there are three impor- 
tant characteristic frequencies that determine the shape 
of the biphoton wave function. The first is the coupling 
effective Rabi frequency f2 e , which determines the two- 
resonance spectrum of the nonlinear susceptibility. The 
second is the linewidth 27,= of the two resonances in the 
nonlinear susceptibility. The third is the full- width- at- 
half-maximum (FWHM) phase-matched bandwidth de- 



termined by the sine function, Acj„ 



2tt x 0. 



where r g is the anti-Stokes group delay time. In time 
domain, they correspond to the Rabi time r r — 2-7r/Q e , 
nonlinear coherence time r e = l/(27 e ) and group delay 
time T g . The group delay time can be estimated from 
T g = L/Vg ~ (27i 3 /|fi c | 2 )0 J D with the optical depth de- 
fined as OD = Nai^L. The competition between r r , 
r e and T g will determine which effect plays a dominant 
role in governing the feature of the two-photon corre- 
lation. Therefore, in the following we will discuss the 
two-photon joint-detection measurement in two regimes, 
damped Rabi oscillation and group delay. 

4. Damped Rabi Oscillation Regime 

Let us first look at the damped Rabi oscillation regime, 
where optical properties of the two-photon amplitude 
(8) are mainly determined by the nonlinear coupling co- 
efficient. This regime requires that the effective cou- 
pling Rabi frequency f2 e and lincwidths 7 e be smaller 
than the phase-matching bandwidth Ao> g or r r > T g 
and l/(27 e ) > r s , so that we can treat the longitudi- 
nal detuning function as $(w) c± 1 and $(r) ~ S(t). 
The biphoton spectral generation rate is proportional to 
\kL\ 2 oc |x^ 3 )L| 2 . Hence both biphoton spectrum in- 
tensity and emission rate are proportional to OD 2 . If 
l/(27 e ) > 2ir/£l e , the two-photon temporal correlation 
clearly exhibits a damped Rabi oscillation with oscil- 
lation period 2ir/£l e ; if l/(27 e ) < 27r/£l e , the second 
and following oscillations are suppressed because of the 
short dephasing time, as experimentally demonstrated in 
[31, 32]. 

Let us consider the case of ft c > I713 — 712 1 first, 
which implies a real effective coupling Rabi frequency 
Q e . The two-photon wave function now is determined by 
Fourier transform of the nonlinear coupling coefficient 
(17), which is 



k(t) = Be 



7eT, 



^H>> 



(24) 



and B = -i T llotTl °T' ■ The heaviside step 
function here is defined as O(r) = 1 for r > 0, and 



0(r) = for r < 0. Combining with Eq. (12), the bipho- 
ton wavepacket becomes 



z e -y e rr e -i(w as +Cl e /2)t aSe -i(w s -a e /2)t 



tf (t oa = t s +r, t.) ~ Lk(T)e-^+^ 
BL 

~Y 

_ e -t(roa»-Oe/2)t o , e -t(w 8 +fie/2)t s -|0/ r \ ^5) 

where vj as + vj s = lu p + lu c and zu as = W31 have been 
used. The physics of Eq. (25) is understood as follows. 
Because the two-photon state is entangled, it cannot be 
factorized into a function of t as times a function of t s . 
\^\ 2 , depends only on the relative time delay r, which 
implies that the pair is randomly generated at any time. 
If the wave packets of the pump and the coupling beams 
arc not plane waves, the term e -H w c+ w p)*£> s would become 
a wave packet with a coherence length determined by the 
pump and coupling. Equation (25) shows also destructive 
interference between the two biphoton channels. The first 
term in the bracket on the RHS of Eq. (25) represents the 
two-photon amplitude of paired anti-Stokes at zu as +£l e /2 
and Stokes at vj s — r2 e /2; while the second term is the 
two-photon amplitude of paired anti-Stokes at Tu as — Q e /2 
and Stokes at w s + O e /2. 

To further see the interference, let us look at the two- 
photon Glauber correlation function, 



G^(T) = ^\BL\ 2 e- 2 ^[l 



3(fi e r)]e(T), (26) 



The physics of Eq. (26) shows the beating (or interfer- 
ence) between two types of paired photons generated 
from the two FWM processes. The two-photon corre- 
lation (26) exhibits damped Rabi oscillations of period 
2n/Q. e . The damping rate is determined by the resonant 
linewidth 27 e in the doublet. At r = 0, G^ vanishes and 
asu 00, G^ also approaches to zero. This indicates 
photon anti-bunching-like effect [G^(0) < G^ 2 \t)]. In 
the counter-propagation geometry, if k c + k p = 0, G^ 2 '(t) 
becomes a symmetric distribution about r = 0, because 
spontaneously emitted paired photons can trigger both 
detectors Di and D2 as experimentally verified in [32, 35]. 
Figure 2(a) shows the first experimental demonstra- 
tion of this type of biphoton source by Balic et al. in 
a cold 87 Rb gas [31]. The solid lines in Fig. 2(al) and 
(a2) are theoretical curves obtained from the Hciscnbcrg- 
Langevin theory. In Fig. 2(bl) and (b2), the curves 
are theoretical simulations using Eq. (26) with the same 
experimental parameters in (al) and (a2). Comparing 
with the Heisenberg-Langevin theory, the state vector 
picture also gives an agreement with the experimental 
data. From Eq. (26), Fig. 2(al) and (bl) correspond to 
the case of r e > r r where the oscillations can be clearly 
observed before damping out; while Fig. 2(a2) and (b2) 
correspond to the case of r e < t t where only the first 
oscillation is observable and other disappear due to the 
fast dephasing rate 7 e = 2i2±2ia. Although in [31] the 
lack of other oscillations in Fig. 2(a2) were interpreted 
due to the group delay, we believe that both cases shown 
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Fig. 2. (color online) Two-photon correlation function in the damped Rabi oscillation regime, (a) Experimental 
data and theoretical curves in (al) and (a2) using the Heisenberg-Langevin theory from Ref. [31]. (b) Theoretical 
simulations using perturbation theory with same parameters from (al) and (a2) and fitted with a vertical scaling 
factor and background accidental coincidence floor. The missing of other oscillation periods in (a2) and (b2) is due 
to the short dephasing time of about 33 ns. The experimental parameters used here are 713 = 714 = 27r x 3 MHz, 
712 = O.6713, OD =11, A p = -7.5 7 i3, and Q p = O.8713. 



in Fig. 2(al)and (a2) are in the damped Rabi oscillation 
regime and the correlation time is well fitted with r e and 
7>. This can be further verified by reducing the optical 
depth by a factor of 1/2, which is corresponding to a re- 
duction of the group delay by a factor of 1/2. However 
the correlation time widthes in Fig. 2(a2) hold the same. 
If fl c < 1 7i3 — 712 1, the effective coupling Rabi fre- 
quency becomes purely imaginary, il e = ±i[3 e . The 
biphoton wave function becomes 



ip(r) ~ Lk(t) = BLe 



sinh(^_)e(r).(27) 



At the weak coupling limit, (3 e — > 713 — 712 and Eq. (27) 
takes the form 



BL 

~2 



t^( T )-» e -'^<W e -7i2r 



e- 713T )6(r). (28) 



This can be understood as follows. Since the coupling 
beam is very weak, the atomic-level structure is not 



changed due to the input power. Two types of bipho- 
tons are generated with the same central frequency but 
are associated with different linewidths, 712 and 713, re- 
spectively. The interference between these two types of 
biphotons gives a manifested exponential decay. In such 
a limit, the EIT disappears and the medium absorption 
decreases the biphoton generation efficiency. This case 
might be observable in the dilute atomic-gas medium 
with very low optical depth (OD -C 1). 

5. Group Delay Regime 

Suggested by Balic el al. [31, 36] and demonstrated by Du 
et al.[33], the group delay regime is defined as r g > t t 
and the (EIT) slow light effect can be used to dynam- 
ically control the biphoton temporal correlation time. 
The group delay condition is equivalent to fl e > Aw 9 , 
i.e., the biphoton bandwidth is determined by the phase- 
matching spectrum. Therefore, in this section we treat 
the third-order nonlinear susceptibility as a constant 
over the phase-matching spectrum. As a consequence, 



the mechanism of double resonance of biphoton gen- 
eration is washed out by the longer group-delay time. 
The two-photon correlation in this regime looks more 
like that with conventional SPDC photons. The coin- 
cidences tends to be a square-like pattern. The biphoton 
wave function can be approximated as tp( T ) — K(jL<&(t) 
where kq is the on-resonance nonlinear coupling con- 
stant. We find that in this regime, the EIT profile mod- 
ifies the biphoton wave function through transparency 
window and slow lig ht. When the EIT bandwidth 
(Auj tr — |ry 2 /(27 13 v / OiD)) is larger than the phase- 
matching bandwidth, the anti-Stokes loss can be ignored 
and the two-photon wave packet approaches a rectangu- 
lar shape. However, when the EIT loss is significant, the 
biphoton wave packet follows an exponential decay. 

When AuJtr > Aw 9 , we can ignore the loss in Eq. (23) 
and rewrite it as 



(2) 



$ H ^ sinc (^)exp( l ^ 



(29) 



'a' ^ r 9' 

Its Fourier transform gives the biphoton wave function 

V>(r) ~ k L$(t) = K V g Ii{T; 0, L/V g ) e - m ^ T . (30) 

The rectangular function II, ranging from r = to 
L/Vg, shows that the anti-Stokes photon is always de- 
layed with respect to its paired Stokes photon because of 
the slow light effect. It is clear that the two-photon corre- 
lation time is determined by the group delay r g = L/V g . 
The photon pair generation rate can be calculated from 
Eq. (18) as 



R=\k \ 



V g L. 



(31) 



Thus even though the spectral generation rate, KoL<J>(u>), 
scales as OD 2 , the total rate of paired counts scales lin- 
early as OD because the bandwidth reduces linearly with 
optical depth [31]. As explained by Rubin et al. [11], 
the rectangular waveform may be understood as follows: 
the paired photons are always produced from the same, 
but unknown, space point in the nonlinear medium. If 
emitted from the front surface, the anti-Stokes photon 
has no delay, and both photons arrive at detectors si- 
multaneously. However, if emitted from the back sur- 
face, the anti-Stokes photon is delayed relative to the 
Stokes photon by r 9 . For conventional SPDC photons, 
the rectangular-shaped biphoton wave packet has been 
well known and its sub-ps correlation time has been di- 
rectly measured in [51]. 

Interestingly, we find that to observe the rectangular 
shape, the condition Aui tr > Aw 9 sets a lower bound for 
the optical depth. Using Aw Jr ~ \fl c \ 2 /(2ji 3 \ / OD) and 
Auj g ~ 2tt/t 9 ~ n\n c \ 2 /(j 13 OD), we obtain OD > An 2 . 
It can be shown that such an optical depth is necessary 
to have the EIT delay-bandwidth product greater than 
two [25]. 

Figure 3(a) shows an ideal rectangular-shape correla- 
tion with a group delay of r g — L/V g . In Fig. 3(b) 
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Fig. 3. (Color online) Two-photon correlation function 
in the group delay regime, (a) Rectangular-shape cor- 
relation for an ideal spontaneous parametric down con- 
verter with a relative group delay L/V g and pair emis- 
sion rate R. (b) Experimental data (red) and theoreti- 
cal curve (solid blue line) obtained from the Heisenberg- 
Langevin theory in Ref . [33] . (c) Propagation delay of an 
anti-Stokes pulse (red) with respect to a reference pulse 
(black). Data are taken from Ref. [33]. (d) Theoreti- 
cal simulation using the perturbation theory fitted with 
same experimental parameters in (b). 



the experimental demonstration of a near-rectangular- 
shape correlation is obtained at the optical depth around 
53 [33], where experimental data are denoted by "+". 
The solid line is theoretical fitting using the Heisenberg- 
Langevin theory. As seen from Fig. 3(b), the correlation 
time of about 300 ns is achieved and is consistent with the 
anti-Stokes pulse delay measurement shown in the inset 



8 



(c) , where the smaller curve represents the delayed signal 
and the higher one is the reference pulse. The experimen- 
tal parameters used here are 713 = 714 = 2n x 3 MHz, 
712 = 0.02 7l3 , OD = 53, A p = 48.67 7 i 3 , n p = 1.16713, 
and fi c = 4.2O713, respectively. The EIT transparency 
width is estimated around 3.63 MHz and the phase- 
matching spectral width about 2.93 MHz. The theoret- 
ical simulation using the perturbation theory [Eq. (13)], 
plotted in Fig. 3(c), also fits very well with the reported 
experimental data. It is found that the exponential-decay 
behavior at the tail is due to the finite EIT loss, which 
alters the correlation function away from the ideal rectan- 
gular shape. In the experiment [33], Du et al. have also 
experimentally verified that the total coincidence counts 
are linearly dependent on the optical depth. 

The sharp peak shown in the leading edge of the two- 
photon coincidence counts in Fig. 3(b) is the first ob- 
served Sommerfeld-Brillouin precursor [52] at the bipho- 
ton level [33, 44]. Mathematically, the integration in the 
two-photon amplitude (13) is similar with the classical 
pulse propagation. In the stationary-phase approxima- 
tion, starting from Eq. (13) one can show that the sharp 
peak is resulted from the interference of biphotons which 
are generated out of the EIT opacity window. The de- 
tailed analysis of the precursor generation at the two- 
photon level has been presented in [44] by using the the- 
ory laid out in Sec. 2. Although by solving the coupled- 
operator equations the Heisenberg-Langevin theory pre- 
dicts the correct pattern, the physics behind the precur- 
sor is not easy to be retrieved from this theory. 

Suppose now the EIT loss is not negligible, i.e, e~~ aL <C 
1. In such a case, the longitudinal detuning function (23) 
can be approximated as 



(a) 



<o„. 



□ <"- 

CO. 



r 



(b) 



{£5 





--> 






R 


a 






a 








Fig. 4. (color online) Biphoton generation in a three-level 
system, (a) Right-angle experimental configuration, (b) 
Three-level double-A energy level diagram. 



the square-wave-like pattern manifested by exponential 
decay will be observed. 

In this section, we are going to briefly show that it is 
possible to switch the two-photon coincidence counting 
rate from photon anti-bunching-like effect to bunching- 
like effect, by manipulating the linear optical responses 
to both Stokes and anti-Stokes fields [43]. We consider a 
case that Stokes and anti-Stokes photons are degenerate 
at both spectrum and polarization so that both Stokes 
and anti-Stokes experience slow light propagation. As 
shown in Fig. 4(a) and (b), we consider paired Stokes and 
anti-Stokes photons produced from a three-level double- 
A scheme in a right angle configuration. In the counter- 
propagation geometry, frequency-phase matching allows 
both Stokes and anti-Stokes fields to reach both detec- 
tors. The two generated fields both may experience slow- 
light effect inside the medium due to the linear optical 
responses. The biphoton wave function now is expressed 
as 



$(w) 



iaL 



The two-photon wave function becomes 



tp(r) ~ K V g e 



-ceVgT —iza a 



(32) 



(33) 



The correlation function |\l/(r)| 2 has an exponential de- 
cay time of T g /(2aL). This means that the anti-Stokes 
photon detected at r is generated at the depth of V g T 
such that the amplitude decreases by e~ aVsT . 

6. Two-Photon Interference: From Anti- 

Bunching to Bunching 

In Sec. 4 and Sec. 5, we have discussed the two-photon 
temporal correlation in two regimes. In general, the two- 
photon wave function is a convolution of both the non- 
linear and linear responses as shown in Eq. (15). In both 
regimes, the correlation functions have non-zero values 
only at r > 0. The observed anti-bunching-like effect in 
the damped Rabi oscillation regime is due to the destruc- 
tive interference between two types of biphoton genera- 
tion. In principle, if the detection system is frequency 
sensitive enough to resolve these two types of bipho- 
tons, the damped Rabi oscillations will disappear and 



rl>{r) ~ [*(r) + S(-r)] * n(r; -L/V g ,L/V g ). (34) 

Even though there is no temporal overlap between k(t) 
and k(—t), their convolution with H(r; — L/V g ,L/V g ) 
mixes them together in the time domain. By manip- 
ulating the group delay, the two-photon interference in 
Eq. (34) can be switched from the destructive to the 
constructive. Alternatively, the coincidence counting 
rate is changed from photon anti-bunching-like effect to 
bunching-like effect. This can be achieved by making the 
indistinguishability of their propagation pathways which 
gives rise to the biphoton interference. The interference 
comes from two indistinguishable paths for paired Stokes 
and anti-Stokes photons. In one possible propagation 
pathway, the pair is created at z\ with Stokes photon 
going to the left and anti-Stokes photon to the right. 
In the other propagation pathway, the pair is created at 
Z2 = —Z\ but with Stokes photon going to the right and 
anti-Stokes photon to the left. The interference between 
the two paths then modifies the two-photon correlation 
by further altering the group delay. The details of this 
two-photon interference are presented in [43]. Recently, 
Wen et al. has proposed a scheme to switch the two- 
photon interference from photon anti-bunching-like effect 



to bunching-like effect by manipulating the linear opti- 
cal responses to the paired photons in a strongly driven 
two- level system [49]. 

7. Conclusion and Outlook 

In conclusion, using perturbation theory, we have dis- 
cussed the nature of the two-photon state and its wave- 
form for narrow-band biphoton generation near atomic 
resonance via FWM parametric process using EIT and 
slow light. In general, the biphoton waveform is de- 
termined by a convolution of the nonlinear response 
(determined by the third-order nonlinear susceptibility) 
and the linear response (determined by the longitudi- 
nal detuning function). Therefore, the two-photon wave 
function can be analytically solved in two regimes, the 
damped Rabi oscillation and the group delay. In the 
damped Rabi oscillation regime, where the biphoton 
wave function is mainly determined by the third-order 
nonlinear susceptibility, the two-photon temporal corre- 
lation exhibits a damped Rabi oscillation because of the 
interference between two types of FWMs. In the group 
delay regime, where the two-photon amplitude is mainly 
determined by the phase-matching condition, the cor- 
relation approaches a rectangular shape manifested by 
an exponential decay. The theoretical simulations agree 
very well with recent experimental results [31, 33] and 
the Heisenberg-Langevin theory [36, 37]. Although the 
discusses are based on a four-level double-A EIT system, 
the analysis can be easily extended to three-level and 
two-level systems. 

Before concluding, we note that a number of other in- 
teresting avenues are presently being explored that are 
closely related with this new type of biphoton source. 
By manipulating the linear optical responses to the gen- 
erated fields, Bell inequality can be tested by two-photon 
beatings using biphotons generated from a two-level sys- 
tem [49]. It is interesting to note that the correspon- 
dence principle fails in describing paired-photon genera- 
tion in a two-level system. As mentioned in [35, 39], it 
is not generally valid to simply replace the coupled-field 
Maxwell's equations by coupled field-operator equations 
without changing the coefficients. But the underlying 
physics is not clear at this moment. 

It will be interesting to explore quantum networks with 
these narrow-band paired photons. For example, the fea- 
sibility of an EIT-based quantum memory can be easily 
demonstrated using current biphoton source by looking 
at the two-photon interference, as suggested in [53]. It 
may be interesting to look at the squeezing properties af- 
ter an EIT system if the input squeezed light is initially 
generated from the atomic-gas medium [54] . As proposed 
in [40, 43] experimental demonstration of polarization- 
entangled narrow-band biphoton generation would be 
useful for quantum computation and communication. 

Another problem is that it is not clear whether the 
maximally time-frequency entangled photon pairs can be 
generated in the room-temperature atomic-gas medium. 
Since all the experimental and theoretical work now are 



focused on cold atomic ensembles, it would be interesting 
to know if the Doppler broadening becomes dominant in 
the process, what kind of feature would be observed in 
the two-photon correlation measurement? By optimiz- 
ing the system, can we still observe the damped Rabi 
oscillations? 

Apart from temporal correlation, a possible applica- 
tion of these narrow-band biphotons may involve the 
transverse correlation [38]. It may be worth studying 
whether this new type of source is useful for quantum 
imaging and lithography. 
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